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Abstract 
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1. Introduction 



In the Hilbert space /2(N), consider the operator J whose matrix representation 
with respect to the canonical basis in /2(N) is the semi- infinite Jacobi matrix 
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where 6„ > and g„ G M for n G N. This operator is densely defined in l2(^) and 
J C J* (see Section |2] for details on how J is defined). 

It is well known that J can have either (1, 1) or (0, 0) as its deficiency indices 
[21 Sec. 1.2 Chap. 4], [31], Cor. 2.9]. By our definition (see Section [2]), J is closed, 
so the case (0, 0) corresponds to J = J*, while (1, 1) implies that J is a non-trivial 
restriction of J*. The latter operator is always defined on the maximal domain in 
which the action of the matrix f 1 1.1 1) makes sense [21 Sec. 47]. 

Throughout this work we assume that J has deficiency indices (1,1). Jacobi 
operators of this kind are referred as being in the limit circle case and the moment 
problem associated with the corresponding Jacobi matrix is said to be indeterminate 
[21 [31]. In the limit circle case, all self-adjoint extensions of a Jacobi operator have 
discrete spectrum ^ST\ Thm. 4.11]. The set of all self-adjoint extensions of a Jacobi 
operator can be characterized as a one parameter family of operators (see Section [2|). 

The main results of the present work are Theorem [1] in Section [3] and Theorem [2] 
in Section m In Theorem [1] we show that a Jacobi matrix can be recovered uniquely 
from the spectra of two different self-adjoint extensions of the Jacobi operator J 
corresponding to that matrix. Moreover, these spectra also determine the parame- 
ters that define the self-adjoint extensions of J for which they are the spectra. The 
proof of Theorem [T] is constructive and it gives a method for the unique reconstruc- 
tion. The uniqueness of this reconstruction in a more restricted setting has been 
announced in [13] without proof. In Theorem [2] we give necessary and sufficient 
conditions for two sequences to be the spectra of two self-adjoint extensions of a 
Jacobi operator in the limit circle case. This is a complete characterization of the 
spectral data for the two-spectra inverse problem of a Jacobi operator in the limit 
circle case. 

In two spectra inverse problems, one may reconstruct a certain self-adjoint oper- 
ator from the spectra of two different rank-one self-adjoint perturbations of the oper- 
ator to be reconstructed. This is the case of recovering the potential of a Schrodinger 
differential expression in ^2(0, 00), being regular at the origin and limit point at 00, 
from the spectra of two operators defined by the differential expression with two dif- 
ferent self-adjoint boundary conditions at the origin [H [51191 fT5l fT71[2n IM 1271128] . 
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Necessary and sufficient conditions for this inverse problem are found in [2^. Char- 
acterization of spectral data of a related inverse problem was obtained in [10]. The 
inverse problem consisting in recovering a Jacobi matrix from the spectra of two 
rank-one self-adjoint perturbations, was studied in [H [121 [HI (SQl [321 [35]. A com- 
plete characterization of the spectral data for this two-spectra inverse problem is 
given in [2U] . 

In the formulation of the inverse problem studied in the present work, the aim is 
to recover a symmetric non-self-adjoint operator from the spectra of its self-adjoint 
extensions, as well as the parameters that characterize the self-adjoint extensions. 
There are results for this setting of the two spectra inverse problem, for instance in 
[T4l [23] for Sturm-Liouville operators, and in [13j for Jacobi matrices. 

It is well known that self-adjoint extensions of symmetric operators with defi- 
ciency indices (1, 1) can be treated within the rank-one perturbation theory (cf. [HI 
Sec. 1.1-1.3] and, in particular, P, Thm. 1.3.3]). Thus, both settings may be re- 
garded as particular cases of a general two-spectra inverse problem. A consideration 
similar to this is behind the treatment of inverse problems in [11] . For Jacobi opera- 
tors, however, the type of rank-one perturbations in the referred formulations of the 
inverse spectral problem are different [B]. Indeed, in the setting studied in [221, one 
has the so-called bounded rank-one perturbations [6l Sec. 1.1]. This means that 
all the family of rank-one perturbations share the same domain. In contrast the 
present work deals with singular rank-one perturbations [6l Sec. 1.3], meaning that 
every element of the family of rank-one perturbations has different domain. Note 
that for differential operators both settings involve a family of singular rank-one 
perturbations. 

The paper is organized as follows. In Section [21 we introduce Jacobi operators, 
in particular the class whose corresponding Jacobi matrix is in the limit circle case. 
Here we also present some preliminary results and lay down some notation used 
throughout the text. Section [3] contains the uniqueness result on the determination 
of a Jacobi matrix by the spectra of two self-adjoint extensions. The proof of this 
assertion yields a reconstruction algorithm. Finally in Section [H we give a complete 
characterization of the spectral data for the two spectra inverse problem studied 
here. 



2. Preliminaries 

Let lfin{N) be the linear space of sequences with a finite number of non-zero 
elements. In the Hilbert space hi^), consider the operator J defined for every 
/ = {fk}'kLi in ^/m(N) by means of the recurrence relation 

iJf)k:=bk-ifk-i + qkfk + bkfk+i, A;GN\{1}, (2.1) 
(^/)i:=gi/i + &i/2, (2.2) 
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where, for n G N, 6ri is positive and g„ is real. Clearly, J is symmetric since it 
is densely defined and Hermitian due to (12. ip and (12. 2p . Thus J is closable and 
henceforth we shall consider the closure of J and denote it by the same letter. 

We have defined the operator J so that the semi-infinite Jacobi matrix (11.11) is 
its matrix representation with respect to the canonical basis {cnj^i in hi^) (see [3l 
Sec. 47] for the definition of the matrix representation of an unbounded symmetric 
operator). Indeed, J is the minimal closed symmetric operator satisfying 



( Jp p A - ( Tp , p ] - nfci\l,/cfci\l\|i| 



We shall refer to J as the Jacobi operator and to (II. ip as its associated matrix. 

The spectral analysis of J may be carried out by studying the following second 
order difference system 

bn-lfn-l + qnfn + bnfn+l = Cfn, n>l, (GC, (2.3) 

with the "boundary condition" 

qifi + hf2 = (fi. (2.4) 

If one sets /i = 1, then /2 is completely determined by (12.41) . Having /i and 
/2, equation (12. 3p gives all the other elements of a sequence {/nj^^Li that formally 
satisfies (12.31) and (12. 4p . Clearly, is a polynomial of ( of degree n — 1, so we 
denote /„ =: P„_i(C). The polynomials Pn{C)y n = 0,1,2, ... , are referred to as the 
polynomials of the first kind associated with the matrix (11.11) [21 Sec. 2.1 Chap. 1]. 
The sequence -P(C) := {A,— i(C)}fcLi is not in //in(N), but it may happen that 

oo 

Y,\Pk{C)f <oo, (2.5) 

A:=0 

in which case P{() G Ker(J* — (I). 

The polynomials of the second kind Q{() ■= {Qfc-i(C)}fcli associated with the 
matrix (II. ip are defined as the solutions of 

bn-lfn-l + qnfn + &n/n+l = C/n , G N \ {1} , 

under the assumption that /i = and /2 = b^^. Then 

Qn-i{C):=fn, VnGN. 

Qn{C) is a polynomial of degree n — 1. 

As pointed out in the introduction, J has either deficiency indices (1, 1) or (0, 0) 
[2l Sec. 1.2 Chap. 4] and [31], Cor. 2.9]. These cases correspond to the limit circle 
and limit point case, respectively. In terms of the polynomials of the first kind, 
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J has deficiency indices (0, 0) if for one G C \ R the series in (12.51) diverges. In 
the hmit circle case (12.51) holds for every C G C [2, Thm. 1.3.2], [311 Thm. 3] and, 
therefore, P(C) is always in Ker(J* — (^/). Another peculiarity of the limit circle case 
is that every self-adjoint extension of J has purely discrete spectrum [3T1 Thm. 4.11]. 
Moreover, the resolvent of every self-adjoint extension is a Hilbert-Schmidt operator 
[231 Lem.2.19]. 

In what follows we always consider J to have deficiency indices (1,1). The 
behavior of the polynomials of the first kind determines this class. There are various 
criteria for establishing whether a Jacobi operator is symmetric but non-self-adjoint. 
These criteria may be given in terms of the moments associated with the matrix, 
for instance the criterion [3ll Prop. 1.7] due to Krein. A criterion in terms of the 
matrix entries is the following result which belongs to Berezans'kii [2, Chap. 1], [TJ 
Thm. 1.5 Chap. 7]. 

Proposition 1. Suppose that sup„gpj < oo and that 'Yl'^=i F~ ^ there is 

G N such that forn>N 

bn-iK+i < bl , 

then the Jacobi operator whose associated matrix is ( fi.il) is in the limit circle case. 

Jacobi operators in the limit circle case may be used to model physical processes. 
For instance Krein's mechanical interpretation of Stieltjes continued fractions [22], 
in which one has a string carrying point masses with a certain distribution along 
the string, is modeled by an eigenvalue equation of a Jacobi operator ^ Appendix] . 
There are criteria in terms of the point masses and their distribution [2], Thm. 0.4 
Thm. 0.5 Appendix] for the corresponding Jacobi operator to be in the limit circle 
case. 

In this work, all self-adjoint extensions of J are assumed to be restrictions of J*. 
When dealing with all self-adjoint extensions of J, including those which imply an 
extension of the original Hilbert space, the self-adjoint restrictions of J* are called 
von Neumann self-adjoint extensions of J (cf. [3, Appendix I], [311 Sec. 6]). 

There is also a well known result for J in the limit circle case, namely, that J is 
simple [21 Thm. 4.2.4]. In its turn this imply that the eigenvalues of any self-adjoint 
extension of J have multiplicity one [3l Thm.3 Sec. 81]. 

Let us now introduce a convenient way of parametrizing the self-adjoint exten- 
sions of J in the symmetric non-self-adjoint case. We first define the Wronskian 
associated with J for any pair of sequences ip = {(pk}'^=i and ip = {ilJk}kLi in h(J^) 
as follows 

iyfc(V3,^) := 6fc(V2fc^/'fc+i - ^fcV^fc+i) , keN. 
Now, consider the sequences f (r) = {ffe(r)}^^ such that, for k E N, 

VkiT):=Pk-iiO) + TQk-iiO), tGM, (2.6) 

and 

VkM := Qk-i{0) . (2.7) 
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All the self-adjoint extensions J(r) of the symmetric non-self-adjoint operator J 
are restrictions of J* to the set [33l Lem. 2.20] 



^^r := {/ = {/fc}r=i e Dom(J*) : lim iy„(t;(r), /) = O} , rGMU{oo}. (2.8) 

n^oo 

Different values of r imply different self-adjoint extensions, so J(r) is a self-adjoint 
extension of J uniquely determined by r [33l Lem. 2.20]. Observe that the domains 
T^r are defined by a boundary condition at infinity given by r. We also remark that 
given two sequences and ip in Dom(J*) the following limit always exists [33j, Sec. 
2.6] 

lim W„(<^,V') =: Woo(<^,V')- 

?i— >oo 

It follows from [311 Thm. 3] that, in the limit circle case, -P(C) and Q{C) are in 
Dom(J*) for every C ^ 

From what has just been said, one can consider the functions (see also [2^ Sec. 2.4 
Chap. 1, Sec. 4.2 Chap. 2]) 

W^{P{Q),P{0)=:D{C), 

W^{Q{<d\P{0)=--B{C). ^ ■ ^ 

The notation for these limits has not been chosen arbitrarily; they are the elements 
of the second row of the Nevanlinna matrix associated with the matrix f 1 1.1 1) and 
they are usually denoted by these letters [21 Sec. 4.2 Chap. 2], [311 Eq. 4.17]. 

It is well known that the functions -D(C) and -B(C) are entire of at most minimal 
type of order one [21 Thm. 2.4.3], [311 Thm. 4.8], that is, for each e > there exist 
constants Ci(e), C2(e) such that 

|/^(C)l<Ci(6)e^l^l, \B{C)\<C,{e)e^\^\. 

If -P(C) is in the following holds 



Thus, the zeros of the function 



D{C) + tB{Q if tGM 
B{C) if T = oo. 



\D{C) + tB{C) if r e M , , 

^r(C) := { ^ ' (2.10) 
\B[Q if r = oo 

constitute the spectrum of the self-adjoint extension J(r) of J. 

A Jacobi matrix of the form fll.ip determines, in a unique way, the sequence 
P{t) = {P„_i(i(:)},^]^, t G M. This sequence is orthonormal in any space L2(M, dp), 
where p is a solution of the moment problem associated with the Jacobi matrix 
(11.11) [2, Sec. 2.1 Chap. 2]. The elements of the sequence {P„_i(t)},^]^ form a basis 
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in L2(M, dp) if p is an N-extremal solution of the moment problem [2, Def. 2.3.3] or, 
in other words, if p can be written as 

p(t) = (E(t)ei,ei), tGM, (2.11) 

where E{t) is the spectral resolution of the identity for some von Neumann self- 
adjoint extension of the Jacobi operator J associated with (11. ip [2, Thm. 2.3.3, 
Thm. 4.1.4]. 

Let p be given by (12. lip , then we can consider the linear isometric operator 
U which maps the canonical basis {e„}J^| in l2{^) into the orthonormal basis 
{Pnit)}n=o in L2{^,dp) as follows 

f/e„ = P„_i, nen. (2.12) 

By linearity, one extends U to the span of {e„}^]^ and by continuity, to all /2(N). 
Clearly, the range of U is all -^2(1^, dp). The Jacobi operator J given by the matrix 
(11.11) is transformed by U into the operator of multiplication by the independent 
variable in L2(M, dp) if J = J*, and into a symmetric restriction of the operator of 
multiplication if J 7^ J* . Following the terminology used in |Tl], we call the operator 
UJU~^ in L2(M, dp) the canonical representation of J. 

By virtue of the discreteness of (j(J(r)) in the limit circle case (here and in the 
sequel, cr{A) stands for the spectrum of operator A), the function pr given by (12. lip , 
with E{t) being the resolution of the identity of J(r), can be written as follows 

Prit) = ^a{\k)'^ , Xkea{J{T)), 

where the positive constant a(Afc) is the so-called normalizing constant of J(r) cor- 
responding to Afc. In the limit circle case it is easy to obtain the following formula 
for the normalizing constants p] Sec. 4.1 Chap 3], [311 Thm. 4.11] 

a{\k) = m\k)\\l^n). \kea{J{r)). (2.13) 

Formula (I2.13p . which gives the jump of the spectral function at A^, also holds true 
in the limit point case, when A^ is an eigenvalue of J [71 Thm. 1.17 Chap. 7]. 

It turns out that the spectral function p^ uniquely determines J{t). Indeed, 
there are two ways of recovering the matrix from the spectral function. One method, 
developed in [16] (see also [32]), makes use of the asymptotic behaviour of the Weyl 
m-f unction 



niriO := 




and the Ricatti equation [IH Eq. 2.15], [32l Eq. 2.23], 

6>(")(C) = gn-C-^r4— , neN, (2.14) 



7 



where mf" (C) is the Weyl m-function of the Jacobi operator associated with the 
matrix (11. ip with the first n columns and n rows removed. 

The other method for the reconstruction of the matrix is more straightforward 
(see [7i Sec. 1.5 Chap. 7 and, particularly, Thm. 1.11]). The starting point is the 
sequence {t'^j^Q, t G M. From what we discussed above, all the elements of the se- 
quence {t'^I^Q are in L2(M, dpr) and one can apply, in this Hilbert space, the Gram- 
Schmidt procedure of orthonormalization to the sequence {t^}tLo. One, thus, ob- 
tains a sequence of polynomials {Pfc(t)}^Q normalized and orthogonal in L2(M, dpT-). 
These polynomials satisfy a three term recurrence equation Sec. 1.5 Chap. 7], [311 
Sec. 1] 

tPk-i{t) = hk-iPk^2{t) + qkPk-i{t)^hkPk{t), ken\{l], (2.15) 
tPo(t) = giPo(t) + &iPiW, (2.16) 

where all the coefficients {k G N) turn out to be positive and {k G N) are 
real numbers. The system (12.151) and (12.161) defines a matrix which is the matrix 
representation of J. 

After obtaining the matrix associated with J, if it turns out to be non-self- 
adjoint, one can easily obtain the boundary condition at infinity which defines the 
domain of J{t). The recipe is based on the fact that the spectra of different self- 
adjoint extensions are disjoint [21 Sec. 2.4 Chap. 4]. Take an eigenvalue. A, of J{t), 
i.e., A is a point of discontinuity of Pr or a pole of m,-. Since the corresponding 
eigenvector P(A) = {Pfc-i(A)}^i is in Dom( J(r)), it must be that 

Woo{v{t),P{X)) = 0. 

This implies that either Woo (<5(0), P(A)) = 0, which means that r = oo, or 

Woo(P(0),P(A)) 

Woo(g(o),p(A)) ■ 

Notation We conclude this section with a remark on the notation. The elements 
of the unbounded set (T(J(r)), r G M U oo, may be enumerated in different ways. 
Let cr(J(r)) = {\k}k£K, where is a countable set through which the subscript 
k runs. If cr(J(r)) is either bounded from above or below, one may take K = N. 
If (j(J(r)) is unbounded below and above, one may set K = Z. Of course, other 
choices of K are possible. Since the particular choice of K is not important in 
our formulae, we shall drop K from the notation and simple write {Xk}k- AH our 
formulae will be written so that they are independent of the way the elements of 
a sequence are enumerated, so our convention for denoting sequences should not 
lead to misunderstanding. Similarly, we write J2kyk instead of '^f.^j^Vk, and the 
convergence of the series to a number c means that for any sequence of sets {Kj}j^^, 
with Kj C Kj+i C K, such that IJ^i -^i — sequence {^^.g^^, tends 

to c whenever j ^ oo. 
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3. Unique reconstruction of the matrix 

In this section we show that, given the spectra of two different self-adjoint ex- 
tensions J{ti), J{t2) of the Jacobi operator J in the hmit circle case, one can 
always recover the matrix, being the matrix representation of J with respect to the 
canonical basis in /2(N), and the two parameters ri, T2 that define the self-adjoint 
extensions. It has already been announced [T3| Thm. 1] that, when ri,r2 G M and 
Ti 7^ ^2, the spectra cr(J(ri)) and 0"(J(r2)) uniquely determine the matrix of J and 
the numbers ti and T2. A similar result, but in a more general setting can be found 
in [m Thm. 7]. 

Consider the following expression which follows from the Christoffel-Darboux 
formula % Eq.1.17]: 

n-l 

E^^'(0 = {Pn-liOP'niO " i^n(C)P^l(C)) = W^n (P(C) , ^'(C) ) • 
k=0 

It is easy to verify, taking into account the analogue of the Liouville-Ostrogradskii 
formula [2^ Eq. 1.15], that 

w.(p(C),P'(C)) = Wn{P{o),P{0)Wn{Q{o),P'{0) 
-wMo),PiO)Wn{P{o),p'{0). 

Thus, 

Y.p^(o = wupio,p'{0) 

k=0 

= D{OB'{C)-B{C)D'{0. 

Indeed, due to the uniform convergence of the limits in (12.90 [21 Sec. 4.2 Chap. 2], 
the following is valid 

B\O = W^{Q{0),P\O) 

Now, a straightforward computation yields (ri,r2 G M, ri 7^ T2) 

9^.,(C)9^;(C) - ^;(C)^..(C) = (r2 - n) [DiOB'iO - BiOD'iO] . 
On the other hand one clearly has 

^n(C)^'oo(C) - ^h;(C)^oo(C) = DiOB'iO - BiOD'iO , new. 
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Hence, 



oo 



^n(c)^;,(c)-9^;,(c)^..(c) 



a 



(0 :=E^'(0 



T2 - n 



(3.1) 



fc=0 



^n(C)^'oo(C) -^;(C)^oo(C) ri G M. 



It follows from fl2.13p that the values of the function a{() evaluated at the points of 
the spectrum of some self-adjoint extension of J are the corresponding normalizing 
constants of that extension. 

The analogue of (13.11) with Ti ^ and ri,r2 G M, for the Schrodinger operator 
in L2(0, cxd) being in the limit circle case is [HI Eq. 1.20]. Formula [HI Eq. 1.20] 
plays a central role in proving the unique reconstruction theorem for that operator 
[T^ Thm. 1.1]. The discrete counterpart of [2], Thm. 1.1] is [T3l Thm. 1]. It is 
worth mentioning that the reconstruction technique we present below is also based 



It is well known that the spectra of any two different self-adjoint extensions of J 
are disjoint [2i Sec. 2.4 Chap. 4]. One can easily conclude this from (13.11) . Moreover, 
the following assertion holds true. 

Proposition 2. The eigenvalues of two different self-adjoint extensions of a Jacobi 
operator interlace, that is, there is only one eigenvalue of a self-adjoint extension 
between two eigenvalues of any other self-adjoint extension. 

Remark 1. One may arrive at this assertion via rank-one perturbation theory, in 
particular by recurring to the Aronzajn-Krein formula [30l Eq. 1.13]. Nonetheless, 
we provide below a simple proof to illustrate the use of (13.11) . The proof of this 
statement for regular simple symmetric operators can be found in [T8l Prop. 3.4 
Chap. 1]. 

Proof. The proof of this assertion follows from the expression (13. ip . It is similar to 
the proof of [21, Thm. 1.2.2]. 

Note that (I2.10p implies that the entire function 9^t(C)) ^ MU{oo}, is real, 
i.e., it takes real values when evaluated on the real line. Let Xk < Xk+i be two 
neighboring eigenvalues of the self-adjoint extension J{t2) of J, with r2 G MU {oo}. 
So Xk, Xk+i are zeros of 91^-2 and by (13. ip these zeros are simple. Since yi'^^{Xk) 
and yi'^^{Xk+i) have different signs, it follows from (13. ip that 9^ri(Afc) and 9^T-i(Afc+i) 
(ri G MU {oo}, Ti 7^ T2) have also opposite signs. From the continuity of on the 
interval [A^, A^+i], there is at least one zero of in {Xk, A^+i). Now, suppose that 
in this interval there is more than one zero of IHt-i , so one can take two neighboring 
zeros of 9^,-^ in {Xk,Xk+i). By reproducing the argumentation above with ti and 
T2 interchanged, one obtains that there is at least one zero of ^K^-a somewhere in 
{Xk, Afc+i). This contradicts the assumption that A^ and Xk+i are neighbors. 



The assertion of the following proposition is a well established fact (see, for in- 
stance [231 Thm. 1]). We, nevertheless, provide the proof for the reader's convenience 



on 



□ 



10 



and because we introduce in it notation for later use. Note that a non-constant en- 
tire function of at most minimal type of order one must have zeros, otherwise, by 
Weierstrass theorem on the representation of entire functions by infinite products 
[25l Thm. 3 Chap. 1], it would be a function of at least normal type. 

Before stating the proposition we remind the definition of convergence exponent 
of a sequence of complex numbers (see Sec. 4 Chap. 1]). The convergence ex- 
ponent pi of a sequence {i'k}k of non-zero complex numbers accumulating only at 
infinity is given by 



pi := inf < 7 G M : lim V . — < oo > . (3.2) 

I r->oo ' Wk\ ' 

\vk\<'r 

We also remark that, as it is customary, whenever we say that an infinite product 
is convergent we mean that at most a finite number of factors may be zero and 
the partial product formed by the non- vanishing factors tends to a number different 
from zero p^, Sec. 2.2 Chap. 5]. 

Proposition 3. Let f{() be an entire function of at most minimal type of order 
one with an infinite number of zeros. Let the elements of the sequence {i'k}k, which 
accumulate only at infinity, be the non-zero roots of f , where {z/fcjfc contains as many 
elements for each zero as its multiplicity. Assume that m G N U {0} is the order of 
the zero of f at the origin. Then there exists a complex constant C such that 

c 



/(C) = Cr lim n 1 - - ' (3-3) 

where the limit converges uniformly on compacts of C 

Proof. The convergence exponent pi of the zeros of an arbitrary entire function does 
not exceed its order [251 Thm. 6 Chap. 1]. Then, for a function of at most minimal 
type of order one, pi < 1. According to Hadamard's theorem [25l Thm. 13 Chap. 1], 
the expansion of / in an infinite product has either the form: 



/(C) = Ce^^' lim n ^ ' ^' ^ ^ 



(3.4) 



Wk\<r 

if the limit 



lim V (3.5) 

r^oo ^ \l>k\ 
Wk\<r 



converges, or 



/(C) = C"^e<+'^ lim 1\g(^;i] , c,deC 



(3.6) 

kfcl<''" 

if (13.51) diverges. We have used here the Weierstrass primary factors G (for details 
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see [25] Sec. 3 Chap. 1]). Let us suppose that the order is one and (13.51) diverges, 
then, in view of the fact that / is of minimal type, by a theorem due to Lindelof 
Thm. 15 a Chap. 1], we have in particular that 



lim y Uj^^ = -c. 

kfel<i" 



This implies the uniform convergence of the series hnir-^oo XI 1 1/^.1 <r ^ compacts 

of C. In its turn, since pi = 1, this yields that lim^^oo n|i/fe|<r (-'- ~ ^) uniformly 
convergent on any compact of C. Therefore, 

lim \\ g( —■ 1 I = e-'^ lim TT ( 1 - ^ ) . 

Wk\<r \i'k\<r 

Thus, (13.61) can be written as (13.31) . 

Suppose now that the hmit (13.51) converges. If the order of the function is less 
than one, then, by [23 Thm. 13 Chap. 1], one may write (13. 4p as (13.31) . If the order 
of the function is one, by [25, Thm. 12, Thm. 15 b Chap. 1], one concludes again 
that (13.41) can be written as (13.31) (cf. Thm 15 in the Russian version of [25] or, 
alternatively, [26l Lect. 5]). 

□ 

Let {A„(r)}„ be the eigenvalues of J(t). In view of the fact that ^r{() is an 
entire function of at most minimal type of order one, by Proposition [3], one can 
always write 



9=1,(0 =aC'Mim n fl-TTr)' reRU{oo}, 



(3.7) 



0<|Afc(T)|<r 



where G R \ {0} and 6r is the Kronecker delta, i. e., 6r = I if r = 0, and 6r = 
otherwise. The limits in (13. 7p converge uniformly on compacts of C Note that when 
r = we have naturally excluded Afc(O) = from the infinite product. 

When writing (13. 7p . we have taken into account, on the one hand, that 9^o(0) = 0, 
which follows from (12.101) and the definition of the function D, and on the other, 
that different self-adjoint extensions have disjoint spectra (see Section [2]). 

Now, let us consider the following expressions derived from the Green's formula 
[21 Eqs. 1.23, 2.28] 



D{O = Cj2P,{0)Pk{O, (3.8) 

fc=0 

00 

B{O = -l + CY,Qk{0)Pk{O- (3.9) 



fc=0 
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Again we verify from (13.81) that -D(O) = 0, while from (13. 9p we have -B(O) = —1. 
Therefore 9^t(0) = — r for every r G M, and 9^oo(0) = — 1. Thus, Cr = — r provided 
that r e M and r ^ 0, and Coo = — 1- 

To simphfy the writing of some of the formulae below, let us introduce -Rr(C) •= 
that is. 



RriC) := C'^ lim J] " YT^) ' r G M U {cx)} 

0<|Afe(r)|<r ^ '^^ 

where 6r is defined as in (13.71) . 

Due to the uniform convergence of the expression 



(3.10) 



d 



n 



0<|Afc(T)|<r 



C 



A, r 



as r — > oo 



one has 



[A,(r)^-Mim n i^-W)) 



0<|Aft(T)|<r 



A,(r) = 



(3.11) 



By (13. ip and (13. lip , one obtains 



Co = a(0). 



(3.12) 



Theorem 1. Let Ti,T2 G MU{oo} with ti ^ T2. The spectra {AA:(ri)}fc, {\k{T2)}k of 
two different self-adjoint extensions J{ti), J{t2) of a Jacobi operator J in the limit 
circle case uniquely determine the matrix associated with J, and the numbers ti and 
r2- 

Proof. For definiteness assume that ti ^ 0, in other words that the sequence 
{^k{Ti)}k does not contain any zero element. 
By (13.11) . we have 



a(A,(ri)) = MRr,i\kin))R'(\kin)) 



(3.13) 



where 



M = ( 



if n,T2eR\{0} 

T2 if Ti = OO , T2y^0 

— Ti if r2 = oo 

[-Co if T2=0 



(3.14) 
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Now, since {a{\k{Ti))}k are the normalizing constants of J{ti) we must have 

1 1 1 



1 



y 1 =J-y 



^ a(Afc(ri)) M ^ i?.,(A,(ri))i?;^(A,(ri)) " 



Therefore 



Thus, M is completely determined by the sequences {Afc(r2)}A: and {\k{ji)}k- In- 
serting the obtained value of M into (13.13^ one obtains the normalizing constants. 
Having the normalizing constants allows us to construct the spectral measure for 
J(ri). Then, by standard methods (see Section [2]), one reconstructs the matrix as- 
sociated with J and the boundary condition at infinity ti. From the value of M and 
Ti one obtains r2, by using the first three cases in fl3.14p . When G {\k{.T2)}k^ one 
does not use ( 13.14p . since it is already known that r2 = 0. 

□ 

Remark 2. Note that the proof of Theorem [1] gives a reconstruction method of the 
Jacobi matrix. Although mentioned earlier, we also remark here that the assertion of 
Theorem [H for the case of Ti, r2 G M, was announced without proof in [T3l Thm. 1]. 



4. Necessary and sufficient conditions 

In this section we give a complete characterization of our two-spectra inverse 
problem. We remind the reader about the remark on the notation at the end of 
Section [2l 

First we prove the following simple proposition related to the converse of Propo- 
sition [31 

Proposition 4. Let {vk}k be an infinite sequence of non-vanishing complex numbers 
accumulating only at oo, and whose convergence exponent pi does not exceed one. 
Suppose that the infinite product 




converges uniformly on any compact of C. Then this product is an entire function 
of at most minimal type of order one if either 5\) converges or if 5\) diverges 
but the following holds 

lini!^ = 0, (4.2) 

r^oo f 

where n{r) is the number of elements of {i'k}k in the circle \(\ < r. 

Proof. Clearly, by the conditions of the theorem, one can express (14.11) in terms of 
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canonical products ^E, Sec. 3 Chap. 1] either in the form 

lim TT (1-^1 = lim TT g(—;0] (4.3) 



r— ►oo 

Wk\<r ^ kfe|<'" 



whenever ( 13.51) converges, or in the form 

lim TT ( 1 - ^ ) = e'^^ lim TT G ( —; 1 ) (4.4) 



r^oo 



otherwise, where 

lim Yl ^k^ = -^- (4-5) 

Wk\<r 

In the case (14.31) . in which the genus of the product is less than the convergence 
exponent of {i/^lfc, it is clear that (14. ip does not grow faster than an entire function of 
minimal type of order one. Indeed, by Thm. 7 Chap. 1], the order of a canonical 
product is equal to the convergence exponent, so when pi < 1 the assertion is 
obvious. For pi = 1 the statement follows from [251 Thm. 15 b Chap. 1]. 

If we have the representation (14. 4p . then pi = 1. By [251 Thm. 7 Chap. 1], the 
canonical product has order one. Since the product of functions of the same order 
is of that same order, the order of (14. ip is one. Then, the assertion follows from [25| 
Thm. 15a Chap. 1] due to (S2D and fOD . 

□ 

Before passing on to the main results of this section, we establish an auxiliary re- 
sult which is related to part of the proof of Theorem 1 in the Addenda and Problems 
of [1 Chap. 4]. 

Lemma 1. Consider an infinite real sequence and a sequence {otj}j of positive 

numbers such that 

^2m 



< OO for all m = 0,1, . . . 

ni ■ 



a, 
J ■' 



Let ^ be an entire function of at most minimal type of order one whose zeros, {Kj}j, 
are simple, and such that 

|5^(it)|->oo as t^±oo, teR. (4.6) 

// 

^ (1 + «:|)[^?'(«:,)P 

then 

m 
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is absolutely convergent for m = 0, 1, . . . , and the absolutely convergent expansion 



holds true for all ( E C \ {Kj}j- 

Proof. The absolutely convergence of (14.71) follows from 



,,.m 



1 + K'^d'{>^j) 



/a,- 



< 



E 



a,- 



< oo 



Construct the function 



MO := 1-5(0^5^ 



5'(k,)(C-k,) 



where the series is absolutely convergent in compact subsets of C \ because 
of (14.71) . Clearly, h{K,j) = for any j. Moreover, it turns out that h is an entire 
function of at most minimal type of order one. To show this, first consider the case 
when J2j Is P^ < cxD. Here, by what we have discussed in the proof of Proposition [3l 
the function d{C)/{C ~ '^^^ expressed by a canonical product of genus zero. 
It follows from [251 Lem. 3 Chap. 1] (see also the proof of [25i. Thm. 4 Chap. 1]) that, 
on the one hand. 



max 

|C|=r 



(C - f^j 



< exp (C(a)r") , pi < a < 1 



for any r > provided that pi < 1. If, on the other hand, pi = 1, then for any 
e > 0, there exists Rq > such that 



max 

|C|=r 



(C - 



< exp (er) 



(4.8) 



for all r > Rq. Hence, in any case, we have the uniform, with respect to j, asymptotic 
estimation (14.81) when < oo. 

Suppose now that IsT^ ~ ^^^^ case, as was shown in the proof of 
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Proposition [3l 



m 

(C - 



K 



|Kfe|<r 



where 



lim } = -c. 



(4.9) 



|Kfc|<r 



On the basis of the estimates found in the proof of |25l Thm. 15 Chap. 1] (see in 
particular the inequahty next to [251 Eq. 1.43]), one can find independent of j, 
such that 



max 

lC|=r 



(C - 



< exp 



c+ E f^k' 

\Kk\<r 



CMimsup^ + e^ + o(- 

r^oo r / \r 



for all r > Ri and e > (see the definition of n{r) in the statement of Proposition Hj). 
Note that if < r, then the above inequality follows directly from the inequality 
next to [25', Eq. 1.43]. If \hij\ > r, the same inequality holds due to 



< 



+ 



Since ^ does not grow faster than a function of minimal type of order one, by 
Thm. 15 a Chap. 1], one again verifies that, for any e > 0, fl4.8p holds for all r greater 
than a certain R2 depending only on the velocity of convergence in the limits (14. 9 p 
and flT2D . 

Thus, one concludes that, for any e > 0, there is R> such that 



max 

[C|=r 



< 



T — 



, ^ , , , IIlcLX 

\C\=r 



(C - ^j) 



< exp(er) 



for all r > R, which shows that h is an entire function of at most minimal type of 
order one. 

Now, the function h/'S is also an entire function of at most minimal type of order 
one [251 Cor. Sec. 9 Chap. 1]. By the hypothesis (14. 6p . 



0. 



which implies that h/^ = (see Corollary of [2S1 Sec. 14 Chap. 1]). 



□ 
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Theorem 2. Let {Xk}k o-nd {^k\k be two infinite sequences of real numbers such 
that 

a) {Afcjfe n {fik}k = 0- For definiteness we assume that ^ {Xk}k 

b) the sequences accumulate only at the point at infinity. 

c) Xk 7^ Xj, Hk 7^ jJ-j for j. 

Then there exist unique ti,T2 G MU {oo}, with ti ^ 0, ti ^ T2, and a unique Jacobi 
operator J ^ J* such that {Xk}k — '7{J{ti)) o-nd {nk}k — c(J(t2)) if and only if the 
following conditions are satisfied. 

1. The convergence exponents of the sequence {Xk}k, and of the non-zero elements 
of do not exceed one. Additionally, if 

lim y —— = 00 , require that lim = , 

^^°°,rrt \Xk\ r 

and if 



lim y - — - = 00 , require that lim — — = , 



Q<\lJ.k\<r 



where n\{r) andn^{r) are the number of elements of {Xk}k o-i^d {^J^k}k, respec- 
tively, in the circle \C, \ < r. 



2. The limits 



XkJ ^^~.^i^.:^^A i^k 



\Xk\<r ^ ' 0<|Mfc|<r 

converge uniformly on compact subsets of C, and they define the functions 

n,(0 lim n - -) ' (4.11) 

where 5 = 1 if & {lJ'k}k, o,nd 5 = otherwise. 

3. All numbers 7l^{Xj)7l'x{Xj) have the same sign for all j. The same is true for 
the numbers TZ\ {/ij ) 7?.^ ( A* j ) • 
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4- For every m = 0,1,2,... the series below are convergent and the following 
equalities hold 

\m . m 

-7? (XATPL (\-\ ~ ~ 2-^ 



^ 7^^(A,)7^l(A,) ^ 7^A(/x,)7^^/i,o 



5. The series 

diverge either to — oo or +cxd. 

6. The series 



V '"^^^^^'^ and V — 

^ (1 + A|)7^',(A,) Z.(i + ^|)7^.(^^.) 

are convergent. 

Proof. We begin by proving that if {Xk}k and {/ifcjfc are, respectively, the spectra of 
the self-adjoint extensions J{ti) and J{t2) of a Jacobi operator J, then conditions 
QH^hold true. 

Since {Xk}k = o"(J(ri)) and {/ifcjfc = o'{J{'T2)), the functions 9^^^ and given 
by (12.101) . have the sequences {Xk}k and {/ifcjfc, respectively, as their sets of zeros. 
These functions do not grow faster than an entire function of minimal type of order 
one. By Proposition [3] (see (13.71) ). the limits 

C\ TT A C 



n 1 - f ' n 1 - 



|Afc|<r ^ ' 0<|/ifc|<r 



converge uniformly on compacts of C. This is condition [H Moreover, by (13. 6p and 
[25l Thm. 15 a Chap. 1], condition [I] holds. 

The functions Tlx and 71^, given by (I4.10p and (14.111) . coincide with R^, given 
by (13.101) . with r = Ti and r = T2, respectively. Thus (I3.13P is rewritten as follows 

a(A,) = M7^^(A,)7^',(A,), (4.12) 

where M is given by (13.141) . Analogously, 

a(/i,) = -M7^A(^,)7^;(/x,). (4.13) 

On the basis of the positiveness of the normalizing constants, from (I4.12p and 
(I4.13p . we obtain condition 

From what we discussed in Section [2] all the moments exist for the spectral 
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functions of J(ri) and J{t2), which are, respectively, 

V-^ and (4.14) 

Hence the series in both sides of condition [7] are convergent for m G N U {0}. 
Moreover, the spectral functions fl4.14p are solutions of the same moment problem 
associated with J (see the paragraph surrounding fl2.1ip ). therefore the equality of 
condition [7] holds. 

Theorem 1 in the Addenda and Problems of [21 Chap. 4] tells us that 

y^^ = +oo (4.15) 

is a necessary condition for the sequences {Xj}j and {a{\j)}j to be the spectrum of 
J(ri) and its corresponding normalizing constants. Thus, substituting (14.121) into 
f l4.15p . one establishes the divergence of the first series in condition O Similarly, 

must hold, which, by fl4.13l) . implies the divergence of the second series in condition 

By the same theorem in [2] mentioned above, and taking into account (14.121) and 
(I4.13p . one obtains the convergence of the series in condition 

Let us now prove that the conditions are sufficient. Using condition [H and 
the convergence of the series in the left hand side of condition [7] with m = 0, we 
define the real constant 

and the sequence of numbers 

a, := Mn^{\^)Tl!^{\j) (4.17) 

By condition and (I4.16p . it follows that aj > for all j. Moreover, (I4.16p and 
fgrrp imply that J^j = 1- 

With the aid of the sequences {\k} and {ofc} define the function p : R — IR+ as 
follows 

Pit) ■= E (4.18) 

Consider the self-adjoint operator of multiplication Ap by the independent variable 
in L2(M, dp). We show below that this operator is the canonical representation (see 
Section [2]) of a self-adjoint extension of a Jacobi matrix in the limit circle case. The 
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proof of this fact is similar to the proof of Theorem 2 in Addenda and Problems of 
[21 Chap. 4]. Note, however, that our conditions are slightly different. 
Consider a function 9k(t) G L2(M.,dp) such that 

9ki\j) = V^kSkj . (4.19) 

Clearly, 6k{t) is the normalized eigenvector of Ap corresponding to A^. Let (f(t) G 
L2(ffi, c/p) be such that 

Taking into account (14.171) . it is clear that the convergence of the first series in 
condition ensures that ip{t) is indeed an element of L2{^,dp). Define 

D:={^e L2(M, dp):^= (A, + i/)" V, ^ e L2(M, dp), ^ ± ^} . (4.21) 

Since D C Dom{Ap), we can consider the restriction of Ap to the linear set D. Let 
us show that this restriction is a symmetric operator with deficiency indices (1,1). 
First we verify that D is dense in L2(M, dp). Suppose that a non-zero r] G 1^2(1^5 dp) 
is orthogonal to D. This would imply that there is a non-zero constant C G C such 
that 

ri = C{Ap~iI)ip. 

Therefore, 



(^7, Oj)L2iR,dp) = 



7^UA,) ' 

whence we easily conclude that 77 G L2(M, rfp) would contradict condition O 

Consider now the restriction of Ap to the set D, denoted henceforth hy Ap f^i, 
and let us find the dimension of ker((74p f^i)* — il) which is characterized as the set 
of all oj G 1^2(1^5 dp) for which the equation 

{{Ap + U)i,uj)L2{W4p) = 

is satisfied for any ^ G -D. It is not difficult to show that any such uj can be written 
as follows 

uj = Cip, ^ C G C. 

Hence dimker((y4p \£,)* — il) = 1. Analogously, it can be shown that the dimension 
of ker((74p f/j)* + il) also equals one. Indeed, if G ker((y4p l^)* + il) then, up to 
a complex constant u = {Ap — iI){Ap + il)~^(p . 

Now we show that Ap I'd is the canonical representation of a Jacobi operator in 
the limit circle case and Ap is the canonical representation of a self-adjoint extension 
of this Jacobi operator. We proceed stepwise. 

I. We orthonormalize the sequence of functions {t"'}'^^^ with respect to the inner 
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product of L2{M.,dp). Note that condition guarantees that all elements of 
the sequence {t"}J^Q are in L2(M, dp). We obtain thus a sequence of poly- 
nomials {Pn-i{t)}'^^i which satisfy the three term recurrence equation fl2.15p 
and (12.161) . where all the coefficients bk {k G N) turn out to be positive and qk 
{k G N) are real numbers. 

II. We verify that the polynomials are dense in L2(M, dp), so the sequence we have 
constructed is a basis in L2(M, dp). Note first that the function 7^a(C) is entire 
of at most minimal type of order one. Indeed, this follows from Proposition HJ 
in view of conditions [I] and [B Now, for any element \ko of the sequence {Xk}k, 
we clearly have 



I^A(it)| > 



t G 



This implies that 7l\ satisfies (14. 6p . Hence the function Tlx and the sequences 
{cij}j satisfy the condition of Lemma [H Thus, we have shown the 
convergence of the series 

A™ 

for all m = 0, 1, 2, ... , and that 



^a(C) ^7^',(A,)(C-A 

Taking into account Definitions 1 and 2 of the Addenda and Problems of [21 
Chap. 4], one obtains from Corollary 2 of [2, Addenda and Problems Chap. 4], 
together with conditions and that {A^} is a canonical sequence of nodes 
and {ci^^} the corresponding sequence of masses for the moment problem given 
by {sm}^=o with 

1 . A™ 

n,{x,)n',{x,) ■ (^-^^^ 

Hence, for this moment problem, p is a canonical solution [21 Def. 3.4.1]. By 
definition, a canonical solution is N-extremal and by [21 Thm. 2.3.3], the poly- 
nomials are dense in L2(M, dp). 

111. We prove that the elements of the basis {P„_i(t)}^]^ are in D. From (I4.22p 
and (14.231) . by Lemma 1 of the Addenda and Problems of [21 Chap. 4], one has 
for m = 0, 1, 2, . . . 
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Then, if S{t) is a polynomial 

Six ■) 

Whence it follows that S E D. 

Now, by ( 12.15P and ( 12.16p . it is straightforward to show that U~^Ap \d U (see 
f l2.12p ) is a Jacobi operator in the limit circle case. 

Denote by J the Jacobi operator U~^Ap \o U . On the basis of what was dis- 
cussed in Section |2] one can find ri G (M U {oo}) \ {0} such that the self-adjoint 
operator of multiplication in L2(M, (ip) is the canonical representation of J{ti). Ti 
cannot be zero since then {Xk} should contain the zero. If ^ {/Ufcjfc, we define 

if Ti = oo 

if n = -M (4.25) 
in all other cases, 

and if G {[ik}k simply assign T2 := 0. 

For the proof to be complete it remains to show that {/ifcjfc are the eigenvalues of 
J{t2). To this end we first show that {iik}k are the eigenvalues of some self-adjoint 
extension of J. Let M. = —A4 and define 

aj := Mnx{Hj)n'p{fij) . 

From condition with m = 0, it follows that aj > for any j and YljO-J'^ = 1, and 
that the function p(t) := X]/xfc<t^fe^ ^ solution of the moment problem {s/c}^q 
with Sk given by fl4.24p . Moreover, taking into account conditions [I] and 0, one 
easily verifies as before that the sequences {fJ,j}j and {aj}j, and the function 7^^ 
satisfy the conditions of Lemma [H Therefore, by Definitions 1, 2 and Corollary 
2 of the Addenda and Problems of [H Chap. 4], as well as conditions and it 
turns out that the sequence is a canonical sequence of nodes and {aj}j the 

corresponding sequence of masses for the moment problem given by {s^j^Q with Sk 
satisfying (14.241) . Hence p is a canonical solution of this moment problem. Denote 
by J{t2) the self-adjoint extension of J having p as its spectral function. 

Let us consider now the functions and Rr^ corresponding to J{t2) and J(t^), 
respectively (see (13.101) ). It is straightforward to verify that 

Rr.i>^,) = mr\x,) = ^^'^^'^ ' ^^-^^^ 

where the first equality follows from (13.131) , while the second follows from (13.111) and 
dHZD. By (I3A4D . (I3T5|) . and (g^SD, one easily concludes from flCTD that = ts- 

□ 



r2 := < 00 

\ Tl+M 
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Remark 3. When G {/ifejfc, the signs of the real numbers 7l^{Xj)7l\{Xj) and 
7l\{fj,j)7l'^{fij) are known. Thus, we can write condition as follows 

n^{x,)n'^{x,) < 7^A(/i,)7^;(/i,) > o for aii j . 

This is a consequence of (13.121) and f l3.14p by which we know that in equation f l4.12p 
M = -a(0) < 0. 

Remark 4. Note that, by Proposition [2], conditions CHS imply the interlacing of 
the sequences {Xk}k and {fik}k- 
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